Autoionization of spin-polarized metastable helium in tight anisotropic harmonic traps 
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Spin-dipole mediated interactions between tightly confined metastable helium atoms couple the 
spin-polarized quintet 5 £+ state to the singlet 1 E+ state from which autoionization is highly prob- 
able, resulting in finite lifetimes for the trap eigenstates. We extend our earlier study on spherically 
symmetric harmonic traps to the interesting cases of axially symmetric anisotropic harmonic traps 
and report results for the lowest 10 states in "cigar-like" and "pancake-like" traps with average 
frequencies of 100 kHz and 1 MHz. We find that there is a significant suppression of ionization, and 
subsequent increase in lifetimes, at trap aspect ratios A = p/q, where p and q are integers, for those 
states that are degenerate in the absence of collisions or spin-dipole interactions. 

PACS numbers: 32.70.Jz, 32.80.Pj, 32.80.Dz, 34.20.Cf 
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I. INTRODUCTION 

There is significant interest in the study and control 
of quantum processes involving tightly trapped ultracold 
neutral atoms where the trapping environment modifies 
the collision properties of the atoms P, 0, Q ■ Optical lat- 
tices, with typical trapping frequencies of 10 4 to 10 6 Hz, 
have recently been used to study quantum phase tran- 
sitions and other phenomena in alkali systems [J] and 
there are now proposals to trap metastable helium and 
neon in optical lattices [E, H 0| • Trapped metastable rare 
gas atoms offer exciting experimental detection strate- 
gies to study quantum gases as the large internal energy 
can be released through Penning and associative ioniza- 
tion during interatomic collisions and through ejection 
of electrons when the atoms strike a metal surface. The 
metastable atoms are generally spin-polarized in order 
to suppress the autoionization rate and to thereby attain 
large numbers of trapped atoms. 

Spin-dipole interactions between the confined 
metastable helium atoms couple the spin-polarized 
quintet 5 S+ state to the singlet 1 E+ state from which 
Penning and associative ionization are highly probable. 
This spin-dipole induced coupling results in finite 
lifetimes for the trap eigenstates and has been studied 
for spherically symmetric harmonic traps by [8| using a 
second-order perturbative treatment of the interaction. 

In this paper we extend these calculations to the in- 
teresting case of anisotropic harmonic traps where there 
arises the possibility of enhancement and/or suppression 
of the ionization losses due to interference between the 
anisotropic part of the trap potential and the spin-dipole 
interaction which are both proportional to a second order 
spherical harmonic Y2 m {0,4>). In particular we consider 
traps with axial symmetry about the z axis. These traps 
have m = and an asymmetry parameter /3 which ranges 
from (3 = 2 for a " pancake" trap to f3 = — 1 for a " cigar" 
trap. 

We have calculated the lowest 10 eigenstates of lu/2tt — 



100 kHz and 1.0 MHz traps where Q is the mean angular 
frequency of the trap. We find evidence that the ioniza- 
tion is strongly suppressed and the lifetimes of the trap 
states significantly lengthened for particular values of (3, 
especially for the higher states considered. 



II. TWO-ATOM COLLISIONS IN AN 
ANISOTROPIC HARMONIC TRAP 

A. Hamiltonian 

Consider two trapped atoms with masses Mi and Mi, 
position vectors r% and r 2 , and electronic spin opera- 
tors S\ and Si subject to a central interatomic potential 
V el (r), where r = \r\ = |ri — 7*2) , and a spin-dipole in- 
teraction @| 



H sd (r) 



[3(Si-r)(S 2 -r)-Si-S 2 ]. (1) 



Here r = r/r is a unit vector directed along the internu- 
clear axis and 



Me 

Ms 



(2) 



where a is the fine structure constant, clq is the Bohr ra- 
dius, (Me/Ms) is the electron magnetic moment to Bohr 
magneton ratio and Eh = a 2 m e c 2 is the Hartree energy. 
Provided the trapping potential is harmonic, the two- 
atom Hamiltonian is separable into center-of-mass and 
relative motions. The Hamiltonian for the relative mo- 
tion is then 



H 



2M 



V 2 . + V cl (r) + H sd (r) + V tlap (r) (3) 



where M = M\Mi / (Mi + M2) is the reduced mass of the 
system. 
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The potential for a general anisotropic harmonic trap 



B. Perturbation theory 



W) = f (<& + & + (4) 

where lj x ,oj v and oj z are the trapping frequencies in the 
x, y and z directions, can be resolved into isotropic and 
anisotropic components by introducing the spherical ba- 
sis x^, [i = 0, ±1, where 



x ± = +-^=(x ±iy), x = z. 
Noting that x^ = rC*(0, <fi) where 



(5) 



4?T 

21 + 1 



Yi, 



(6) 



is the modified spherical harmonic then, using the prod- 
uct formula 

f_i \mi 



C(l, 1, 2; mi, m 2 ,mi + m 2 ) 



xc; 



rn.i+rri2 \" ' 



(7) 



where C(j\,j2,j; n%i, m 2 ,m) is the Clebsch-Gordan coef- 
ficient, Vtrap(^) can be written as 

V tiap (r) = V™ p (r) + V££°(r, 9, 4>). (8) 

The isotropic component is 

K?a P W = (9) 
and the anisotropic part is 

V£*°{r, 6, 4>) = V™ p (r) { [C 2 2 (0, cf>) + Cl 2 (d, </>)] 

+0Cg(M)}- (10) 
We have introduced the mean square frequency 



1 



and the anisotropy parameters 

1 



(11) 



, ^=—(2^-^-^). (12) 



For axially symmetric traps uj x — u> y and a = 0. The 
parameter /? then varies over the range — 1 < (3 < 2 where 
(3 = —I corresponds to an infinitely long cigar-shaped 
trap and /3 — 2 to an infinitesimally thin pancake-shaped 
trap. 



The spin-dipole interaction is of 0(a 2 ) and can be 
treated as a perturbation. The unperturbed Hamiltonian 
Hq is then that for ultracold collisions in an anisotropic 
harmonic trap. The two colliding atoms with spin quan- 
tum numbers Si and S 2 are in the total spin state 
\(SiS2)SMs) where Ms is the projection of S onto the 
space-fixed z axis. During the collision they form the 
molecular state 2S+1 A where A is the projection of the to- 
tal electronic orbital angular momentum along the molec- 
ular axis. The adiabatic potential of this state will be 
denoted by V£ s (r). The appropriate channel states are 



111 



|$„) = \TSM s )\lm) 



(13) 



where a = {T, S, Ms,l,m} and \lm) = Yi m (9,<f>) are the 
relative motion states. The label T represents the remain- 
ing quantum numbers Si, S 2 , ■ ■ ■ needed to fully specify 
the channel. The unperturbed Hamiltonian Hq is given 
by 

Ho = + Vf s (r) + V/ trap (r) = + Vg£°(r) 

(14) 

and the eigenstates |\l/„) satisfy 

#o|*n> = E n \V n ) . (15) 

The states |\E' n ) can be expanded in terms of the channel 
basis defined in (fTB")) so that 



(16) 



and we obtain the matrix eigenvalue equation 

^2H a , a R na {r)=E n R na ,(r) (17) 

a 

where, assuming the r-dependence of the channel states 
is negligible, 

H a , a = (*„/|.Ho|$„) 

= h%>5 a ,, a + (* a ,\VSg°\$ a ) (18) 



and 



h*> = (<t a \H s °\<!> a ) 

2M dr 2 
Muj 2 o 



l(l + i)h 2 

2Mr 2 



(19) 



We have distinguished between right eigenvectors |$ a ) 
and left eigenvectors ($ a | , see @ , as Ho is non-Hermitian 
for the singlet state. In the absence of the anisotropy 
term equation l|17p decouples to 



2M dr 2 



1(1 + i)h 2 

2Mr 2 



+ VZ s (r) + 



R a (r) 
E a R a (r) ,(20) 
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where Egj = 2hui(n* + 1/2 + 3/4) are the energy eigenval- 
ues, see [Jf, of an isotropic harmonic trap of frequency Q. 
The quantum number n* = n r — fx replaces the usual 
integer n r > labelling the cigcncncrgics of a three- 
dimensional isotropic oscillator and includes the effect of 
collisions through the quantum defect /i which depends 
on {A,S,l}. 

The channel basis elements in (jTTJ) can be evaluated 
using standard techniques. From (fT0|) we have 

(*-'|V££°|*a) = Sy^{l'm'\V^°\lm) (21) 

where 7 = {T, S, Ms}. Using the Wigner-Eckart theorem 
then gives 



(*.'|VS^°|*.) = VSp(r)P^ p 



(22) 



where 



D Ta = {a[C(l,2J';m,2,m')+C(l,2J';m,-2,m')} 
+0C(l, 2, 1'; m, 0, to')} (l'\\C 2 \\l)S Yn . (23) 

The reduced matrix element is given by 
~ 21 + 1~ 1 



(i'\\c 2 



2V + 1 



C(Z,2,Z';0,0,0). 



(24) 



To second order, the change in energy of the state \^f n ) 
due to the perturbation H sd is 



AE„ = AM 1 ' 



AE£\ 



where 



AElp = (y n \H sd \y n ) 



and 



(gggtjtt m )(tt m |ff sd |tt t 

E n — Em 



(25) 
(26) 

(27) 



The evaluation of the first order correction is straightfor- 
ward. The second order correction is conveniently evalu- 
ated by introducing the operator F n which satisfies 

{9 m \[F n ,H ]\* n ) = (¥ m |fl Bd |¥ B ) (28) 
= {E n -E m )(V m \F n \* n ) (29) 

for m =/= n. Invoking the closure relation 
J2 m l#m)(*m| = I then gives 

AE™ = ($ n \H sd F n \* n ) - (*„|flad|*n)<*»|^|*n>. 

(30) 

Equation ([28]) yields the inhomogeneous equation 

[F n ,H ]\V n ) = {H sd - ABW)|*„) (31) 
so that, expanding F n \^> n ) in the channel basis (fT3|) 

F n \*n) =^2 fna(r)\$ a ) (32) 



gives the coupled equations 

Y}E n 5 a >,a - H a , a ]f na (r) = ^[($ a ,|ff sd |$ a ) 

a a 

-AE^8 a ,. a ]R na {r). (33) 

To evaluate the elements for the spin-dipole interaction 
we note that the interaction may be written as the scalar 
product of two second-rank irreducible tensors, see [g], 

H sd = V sd (r)T 2 -C 2 , (34) 

where T 2 is 

T 2 ee [Si X S\\l = C(l, 1, 2;/j,,i> — M , u) S^S^, 

and C 2 is the second rank tensor formed from the mod- 
ified spherical harmonics. The radial factor is V sd (r) = 
b/r 3 where b = —^/6k/h 2 . We then get 



($ a ,\H sd \$ a ) = V sd (r)D*J a 



(36) 



where, by invoking the Wigner-Eckart theorem, we find 

El 



Dl'a - $M s ,+m>,M s +m (-l) Ms ' Ms 

x C(S, 2, S'; M s , M s , - M s , M s , ) 
xC(l, 2, l'\ m, m' — m, to') 

x(r's'||r 2 ||rs)(z'||c 2 ||0. 

The reduced matrix element of T 2 is given by 
(r'S'||T 2 ||rS) = Sr',rS s{ , Sl 6 S , S2 fi 2 

y.^S 1 {S 1 + l)S 2 {S 2 + l) 

S\ S2 s 
1 1 2 

$1 $2 $ f 



(37) 



(38) 



where the angular momentum coefficient in (|38[) is related 
to the Wigner 9 — j coefficient 



3i 32 3 
fci hi k 

3i 3 2 3 



ee [(2# + l)(2j£ + l)(2j + l)(2fc+l)]* 



(39) 



In terms of the channel basis expansions, the energy 
shifts are given by 




A£«=]r D s * a (R na ,\V sd \R n 

a' ,a 



(40) 



and 



AE& = J2 D S a l a (Rna>\V sd \fna) ~ A£« ^(^al/r, 



(41) 
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III. APPLICATION TO ULTRACOLD 
SPIN-POLARIZED METASTABLE HELIUM 
ATOMS 

We now consider application of the theory to spin- 
polarized metastable helium atoms tightly confined in 
axially symmetric harmonic traps for which the asym- 
metry parameter a in (|f 2[) vanishes. The colliding atoms 
are in the \(SiS 2 )SM s ) = | (11)22) 5 £+ molecular state 
but, noting that the reduced matrix elements for the spin- 
dipole interaction are [ll[ (r'2||T 2 ||r2) = fr'.r h 2 ->Jfj3 
and (r'0||T 2 ||r2) = -8 r >,r tf-s/SjZ ((r'l||T 2 ||r2) van- 
ishes by symmetry of the 9 — j symbol) , then the spin- 
dipole interaction couples this state to the | (11)00) 1 E+ 
state from which Penning and associative ionization pro- 
cesses are highly probable. We model these ionization 
processes from the singlet state using a complex optical 
potential of the form V{r) — iT{r)/2 where T(r) is the 
ionization width. This coupling between the initial quin- 
tet state and the singlet state produces complex energy 
shifts AE n — AE^ C — ij n /2 due to the complex form of 
the singlet potential. The 1/e lifetimes of the trap states 
will be T n = h/j n [l3|. 

The eigenvalue equations (fl~7|) for R na {r) and the cou- 
pled differential equations (I33|) for the perturbed func- 
tions / no (r) were solved using a discrete variable repre- 
sentation and a scaled radial coordinate grid p = r/£ = 
(,t a where £ = yjTi/Mu) is the effective range of the oscil- 
lator ground state and ( and a are scaling parameters (see 
3 for details). We note that the bosonic symmetry of 
the identical metastable helium atoms restricts the scat- 
tering to even partial waves in the expansions (|16[) and 

As input to the problem we require the Born- 
Oppenheimer potentials for both the 5 S+ and X E+ elec- 
tronic states of the metastable helium dimer. For the 5 E+ 
potential we have used the full potential Vp j (r) of Przy- 
bytek and Jeziorski [l2j which includes adiabatic and rel- 
ativistic corrections, adjusted to match the experimental 
binding energy of the least bound (v = 14) state @ . For 
the singlet potential we use a potential constructed 
from the short-range potential of Miiller et al [l3| expo- 
nentially damped onto the quintet potential at long range 
[Til ]. Also required is the ionization width T(r). Two 
forms were used, TmM obtained from a least squares 
fit to the tabulated results in [l3[ and the simpler form 
r G MsM = 0.3exp(-r/1.086) advocated in 



IV. RESULTS AND DISCUSSION 

We have calculated the lifetimes of the lowest 10 trap 
states for axially symmetric anisotropic harmonic traps 
with average trapping frequencies of 100 kHz and 1 MHz. 
We label these states v = 0, . . . , 9. 

With the power law scaling a = 10, three to four digit 
convergence with respect to the radial grid was obtained 



using 500 (450) grid points for w/2tt = 100 kHz (1 MHz) 
and an outer boundary set at 15£. The position of the 
outer boundary was determined by the requirement that 
the states up to at least v — 9 were insensitive to its 
position. A similar three to four digit convergence with 
respect to the number of partial waves included was gen- 
erally obtained with six partial waves (I = 0, 2, ... , 10) 
for asymmetry values in the range —0.9 < (3 < 1.9, where 
A/3 = 0.01. However, for some of the cases where the life- 
times are strongly enhanced (see below), satisfactory par- 
tial wave convergence was not obtained and the results 
reported for these cases should be regarded as indicative 
only. 

The lifetimes of the six lowest states for lo/2tt = 100 
kHz are shown in Figs Q] and as a function of the trap 
aspect ratio A. For "cigar-like" traps with u> x > lo z the 
aspect ratio is 



(1 + /3) 



(42) 



whereas for "pancake-like" traps with oj z > lu x it is A = 
uj z /uj x . 

The results display a strong sensitivity of most life- 
times to the aspect ratio A, especially around the values 
A = p/q where p and q are integers. Under higher reso- 
lution these peaks show a double peak structure centred 
on A = p/q. In Table 1 results are presented for the 10 
lowest states at the (3 values closest to these aspect ratios 
and it is clear that there is an aspect ratio for nearly ev- 
ery state at which the lifetime is greatly enhanced. The 
behavior near A = 1 ((3 = 0) is also interesting. For 
(3 = the states become pure partial waves and those 
states v = 0, 2, 5, 9 which contain I = develop lifetimes 
which increase slowly for small \/3\ . However the other 
states show strongly enhanced lifetimes at very small val- 
ues of (3 ~ O(10 -3 ) either side of j3 = 0. No results are 
reported for the v = 4, 7, 8 states at A = 1 as for these 
cases our program developed numerical instabilities due 
to ill-conditioned matrices at very small (3 < O(10 -4 ). 
We note that the results given in [8j for the / = states 
of a spherically symmetric trap are reproduced for (3 = 
if we use the potential Vg M (r) of Starck and Meyer [ll| . 

To understand this behavior we note that, in the ab- 
sence of collisions, the energy eigenvalues for an axially 
symmetric harmonic trap 

£ no - co11 = nu> m (n a + l) + 7w y (n y + -) + fw z (n z + ^) (43) 



can be written as 



E, 



no-coll 



y/l + (3{2n z + 1 + 2A{n x + n y + 1)] (44) 



where A = uj x /lu z . For A = p/q this becomes 



^no-col! = ^y_JL_ [(2nz + l)q + 2{nx +ny + 1)pV 



(45) 
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FIG. 1: Lifetimes of the six lowest trap states for "cigar-like" traps with u/27r = 100 kHz as a function of trap aspect ratio 
A = lu x /uj z . 




FIG. 2: Lifetimes of the six lowest trap states for "pancake- like" traps with ui/2n = 100 kHz as a function of trap aspect ratio 
A = uj z /uj x . 
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TABLE I: Lifetimes (in seconds) of the 10 lowest trap states for axially symmetric anisotropic harmonic traps of selected aspect 
ratios A and an average trapping frequency of 100 kHz. The upper part of the table is for "cigar-like" traps where A = uj x /u z , 
the lower part is for "pancake-like" traps where A = uj z /lu x . The numbers in parentheses denote powers of 10. For many of 
the cases of very strongly enhanced lifetimes the results are indicative only as satisfactory partial wave convergence was not 
achieved. 



A 


v = 


v = 1 


v = 2 


v = 3 


v = 4 


v — 5 


v = 6 


v = 7 


v = 8 


v = 9 


1 


9.977 


5.76(-6) 


6.714 


1.35(-5) 




5.399 


2.28(-5) 






4.645 


1.25 


10.14 


23.86 


9.604 


30.55 


22.83 


9.618 


34.19 


29.36 


22.80 


10.04 


1.5 


10.50 


22.44 


10.54 


27.49 


22.42 


134.8 


8.592 


27.44 


93.43 


17.31 


1.75 


10.98 


22.53 


11.91 


24.39 


23.90 


34.66 


11.94 


26.20 


33.67 


24.99 


2 


11.51 


23.04 


689.5 


8.540 


384.8 


14.78 


1.144(3) 


941.1 


7.135 


554.3 


2.25 


12.01 


23.67 


40.32 


11.38 


41.24 


22.32 


53.39 


36.95 


11.28 


55.80 


2.5 


12.57 


24.46 


36.83 


12.70 


37.42 


24.29 


47.06 


34.21 


1.76(3) 


10.32 


2.75 


13.01 


25.09 


36.50 


13.96 


33.76 


25.80 


44.35 


33.79 


63.40 


13.22 


3 


13.56 


25.91 


36.76 


82.08 


11.86 


36.17 


30.92 


36.43 


49.55 


16.60 


3.25 


14.02 


26.60 


37.20 


64.36 


12.95 


54.67 


24.51 


58.80 


31.72 


95.24 


3.5 


14.57 


27.44 


37.86 


52.19 


14.49 


47.73 


26.70 


55.57 


31.84 


73.23 


3.75 


15.27 


28.56 


40.11 


44.31 


18.01 


39.35 


24.84 


59.83 


27.37 


66.31 


1.25 


10.16 


11.08 


18.00 


11.24 


19.14 


25.52 


11.14 


19.36 


26.77 


35.07 


1.5 


10.60 


10.95 


21.28 


10.73 


21.73 


2.2(2) 


8.133 


21.34 


206.9 


8.306 


1.75 


11.22 


11.36 


25.85 


10.62 


26.04 


10.90 


32.11 


25.00 


10.44 


34.55 


2 


11.98 


11.99 


4.0(6) 


8.042 


6.1(6) 


8.017 


4.2(3) 


3.7(2) 


6.451 


2.2(2) 


2.25 


12.80 


12.71 


14.08 


22.30 


13.94 


21.50 


14.89 


21.28 


29.26 


14.34 


2.5 


13.67 


13.50 


14.12 


27.25 


13.67 


25.73 


14.09 


23.91 


96.42 


10.96 


2.75 


14.60 


14.36 


14.71 


33.25 


13.67 


29.99 


13.95 


28.57 


13.88 


34.22 


3 


15.55 


15.23 


15.41 


1.2(2) 


11.23 


65.20 


11.73 


39.89 


12.46 


52.50 


3.25 


16.57 


16.19 


16.21 


17.98 


28.07 


17.24 


26.14 


17.64 


22.06 


18.29 


3.5 


17.67 


17.19 


17.06 


17.65 


34.87 


16.86 


30.27 


17.24 


25.13 


17.06 


3.75 


18.79 


18.28 


18.21 


18.53 


43.46 


16.95 


39.90 


16.94 


33.72 


17.01 


4 


19.87 


19.29 


19.12 


19.23 


1.9(2) 


13.82 


82.23 


14.63 


34.53 


17.59 


4.25 


20.82 


20.23 


20.42 


20.95 


23.07 


36.13 


22.20 


26.31 


20.65 


23.55 


4.5 


21.94 


21.22 


20.82 


20.65 


21.71 


41.38 


20.64 


38.00 


20.14 


31.79 


4.75 


23.27 


22.45 


21.88 


21.56 


22.03 


50.10 


20.45 


43.68 


20.16 


35.71 


5 


24.30 


23.40 


22.68 


22.26 


22.50 


64.85 


19.68 


51.71 


19.62 


2.2(3) 



From (|45|) it is clear that, for a given pair of integers 
(p,q), there will be states (n x ,n y ,n z ) that are degen- 
erate. For example, consider the case (p, q) — (3, 2) 
for which the lowest degenerate states (n x ,n y ,n z ) 
are (1, 1, 0), (0, 0, 3), corresponding to the vibrational 
states v = (5,6) and (1, 1, 1), (0, 0, 4) corresponding 
to v = (8,9). Collisions then break the degeneracies 
of the eigenstates formed from even partial waves. 
Inspection of Figs Q] and [2] and Table 1 shows that it is 
these degenerate oscillator states that have the strongly 
enhanced lifetimes. 

The enhancement of some lifetimes at very small (3 can 
also be traced to underlying degeneracies. As (3 — > 0, the 
energies approach the isotropic oscillator eigenvalues 

£ no - coll (/3 = 0) = ^(4n r + 2Z + 3) (46) 
which have the degeneracies (n r ,l) — (0,2), (1,0) corre- 



sponding to v = 1, 2, (0, 4), (1, 2), (2, 0) to v = 3, 4, 5 and 
(0, 6), (1, 4), (2, 2), (3, 0) to v = 6, ... , 9. Again collisions 
break some of these degeneracies, predominantly by rais- 
ing the energies of the I = states v = 0, 2, 5, 9. 

The results for Hj/2ir = 1 MHz show the same general 
features as those for 100 kHz and are not presented here. 

In conclusion, our calculations show that the ionization 
losses are sensitive to the anisotropic trapping environ- 
ment and that significant suppression of ionization can 
occur at particular trap aspect ratios. Although there 
is some dependance of our results on the choice of ion- 
ization width r(r) (the differences in lifetimes obtained 
using the ionization widths Tm and Tgms are remark- 
ably constant with the Toms results lower by 14.2%) and 
satisfactory partial wave convergence could not obtained 
with the limit of six partial waves for some of the strongly 
enhanced lifetimes, the general features displayed by our 
results are robust. 
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